UNCLASSIFIED 


AD  NUMBER 


AD873241 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  Jun  1970. 
Other  requests  shall  be  referred  to 
Commander,  Air  Weather  Service  [MAC] , 
Attn:  DN,  Scott  AFB,  II  62225 


AUTHORITY 


AFWA  ltr,  19  Apr  2005 


THIS  PAGE  IS  UNCLASSIFIED 


AD  S-73-2^ 


Technical  Report  234 


of 

IONOSPHERIC  TOTAL 
ELECTRON  CONTENT 


By 

ALLAN  C.  RAMSAY 

MAJOR  US  A  F 

HEADQUARTERS  AIR  WEATHER  SERVICE 

p p ro c~ej  foe-  P o  b  l-  C.  & C  L ea5<2-  ,  Qj 

Th  i  g  d  o-c  u  m-ont  i  s-6 i»  Ij  j  o  c  t  t  o  s  p  o  ci  o  I  0  *p-ert  t  o  n  t  r  oJs 
-o-b  d-e-a  c  h-i  r  a  rt-s-  m  ittal  to  foreign  go  v-cr  n  m  0  n  t  s 


i+h- 


or  foreign  nelionals  may  bo  made  only  wi 
-p  r-i  o  r  a  p-p  roval  of  Hq,  Air  Weather  Service 


PUBLISHED  BY 


SB tifb 

‘mm 

4  f\v>&  2005" 


AIR  WEATHER  SERVICE  (MAC) 
UNITED  STATES  AIR  FORCE 

JUNE  1970 


June  1970 


Technical  Report  23^ 


PREFACE 


Ionospheric  total  electron  content  has  a  significant  effect  on  time-delay 
measurements  of  trans ionospheric  radio  signals.  Radio -ranging  techniques  used 
by  radar  and  future  satellite  navigation  systems  must  account  for  the  variable 
time-delay  introduced  by  the  ionosphere. 

This  report  presents  a  method  for  specifying  or  predicting  the  total  elec¬ 
tron  content  of  the  undisturbed  mid-latitude  ionosphere  during  the  maximum 
phase  of  the  solar  cycle.  This  report  is  intended  to  benefit  developers  and 
operators  of  transionospheric  range -measuring  systems;  it  is  operationally 
oriented,  and  is  not  intended  as  a  review  paper  on  global  electron  content  or 
as  a  discourse  on  ionospheric  physics. 

I  gratefully  acknowledge  the  cooperation  and  assistance  of  both  the  Air 
force  Cambridge  Research  Laboratories  (particularly  J,  A.  Klobuchar,  in  the 
development  and  refinement  of  this  electron-content  model)  and  Major  E.  W. 
Friday,  Jr.,  Hq  Air  Weather  Service,  who  undertook  the  effort  of  preparing  the 
Appendix  to  this  report. 
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LIST  OF  SYMBOLS 

Basic  component  of  solar  10.7-cm  radio  flux. 

Basic  component  of  "standard"  total  electron  content. 

Constant  of  proportionality  relating  total  electron  content  to 
the  active -region  component  of  solar  radiation. 

Electron  charge  (4.8  x  10~1Q  esu)/or  electrons. 

Frequency  (Hz  or  mHz). 

Critical  frequency  of  F2  region. 

Daily  value  of  solar  10.7-cm  radio  flux  (10"3a  watts/ma/Hz). 
Mean  value  of  F10  over  x  days. 

Height  above  earth's  surface. 

Factor  relating  TEC  at  any  site  to  TEC  at  standard  site 
(Section  F:  MOD  III). 

Factor  relating  TEC  at  any  site  to  TEC  at  standard  site 
(Section  E:  MOD  II). 

Electron  mass  (9*1  x  10~aa  gm). 

Refractive  index. 

Electron  number  density. 

Maximum  number  density  of  electrons  in  F2  region. 

Radius  of  earth. 

Range  error  (meters). 

Total  Electron  Content  (10**  e/m8). 

Zenith  angle;  angle  between  ray  direction  and  the  zenith  at 
ionospheric  penetration  point. 

Slab  thickness  (km). 

Latitude. 

Noon  solar  zenith  angle. 

Angular  frequency  (rad/sec). 


Referring  to  any  point. 

Rtf erring  to  the  standard  point:  ionospheric  penetration 
between  Hamilton,  Massachusetts,  and  ATS-3*  near  40* N  7‘/v. 
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p  -  Referring  to  peak  value  in  diurnal  curve, 

sp  -  Referring  to  peak  diurnal  value  at  standard  point, 

x  -  Number  of  days  over  which  P  is  meaned. 
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A  MODEL  OP  IONOSPHERIC  TOTAL  ELECTRON  CONTENT 


SECTION  A  —  INTRODUCTION 


Organization 

This  report  is  written  for  application  to  a  specific  and  limited  operation¬ 
al  problem:  that  of  accounting  for  ionospheric  time  delay  in  VHF  and  UHP 
transionospheric  range -measuring  systems.  It  provides  only  the  briefest  back¬ 
ground  on  the  ionosphere;  emphasis  is  given  to  the  development  and  presenta¬ 
tion  of  models  of  ionospheric  total  electron  content.  The  models  have  been 
constructed  primarily  for  operational  utility  rather  than  scientific  elegance. 
Seasonal,  diurnal,  and  solar-radiation-related  variation  in  total  electron 
content  are  modeled;  impulsive  changes  related  to  geomagnetic  disturbances  are 
not  modeled,  but  are  described  qualitatively. 

The  number  of  references  used  in  this  report  has  been  minimized,  and  ref¬ 
erences  have  been  restricted  wherever  possible  to  comraonly-available  sources. 
Readers  should  consult  referenced  material  for  in-depth  discussions  of  physical 
processes  and  limitations  of  measuring  techniques,  and  for  more  extensive 
bibliographies. 

The  report  begins  with  a  brief  review  of  the  source  of  range  errors,  a 
description  of  the  earth's  ionosphere,  and  a  description  of  relevant  portions 
of  the  solar  radiation  spectrum.  Section  D  presents  a  basic  means  of  specify¬ 
ing  or  predicting  peek  daily  values  of  total  electron  content  at  a  "standard" 
mid-latitude  location.  "Standard"  refers  to  a  site  for  which  considerable  data 
were  available  for  specifying  peak  daily  values  of  electron  content  and  its 
diurnal  variation.  The  "standard"  for  work  presented  here  is  the  APCRL 
Sagamore  Hill  Radio  Observatory.  Provision  is  mads  for  the  artificial  con¬ 
struction  of  a  typical  diurnal  curve  of  total  electron  content.  Section  E  ex¬ 
tends  the  "standard"  electron-content  peak  value  and  the  empirical  diurnal 
variation  of  Section  D  to  other  locations  through  the  use  of  ionospheric  clim¬ 
atology.  Section  P  again  applies  ionospheric  climatology  to  the  "standard" 
peak- value  prediction  generated  it.  Section  D,  but  additional  assumptions  are 
made  to  eliminate  problems  in  timing  of  the  diurnal  variation.  The  version 
presented  in  Section  P  is  the  recommended  general-purpose  total-electron- 
content  model.  Section  0  describes  areas  of  investigation  that  would  be  moat 
likely  to  improve  the  validity  and  eebraclveness  of  the  model;  slab  thickness 
variations,  solar  cycle  variations,  and  geomagnetic -disturbance  effects  are 
diecueaed. 
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Range  Measurements 

The  determination  of  transmitter-to-target  range  by  measuring  the  transit 
time  of  radio  signals  is  a  technique  that  will  be  used  by  Air  Force  radar  and 
navigation  systems  for  years  to  come.  Range  is  determined  simply  and  directly 
by  multiplying  the  transit  time  of  an  electromagnetic  signal  by  the  velocity 
of  signal  propagation.  However,  where  the  refractive  index  of  the  propagating 
medium  differs  from  unity,  the  velocity  of  the  signal  is  decreased,  and  in¬ 
ferred  ranges  can  be  in  error. 

For  radio  signals  propagating  within  the  troposphere,  range  errors  (dif¬ 
ferences  between  inferred  and  geometric  ranges)  are  independent  of  the  frequen¬ 
cy  of  the  signal  and  reach  a  maximum  of  about  100  meters  [19)*  For  radio  sig¬ 
nals  which  travel  through  the  ionosphere,  however,  range  errors  whose  magnitude 
is  inversely  proportional  to  the  square  of  the  frequency  can  sometimes  be 
measured  in  kilometers.  The  major  range  errors  introduced  in  transionospheric 
VHP-UHF  signals  can  be  accounted  for  if  the  nature  of  the  propagating  medium 
is  known. 

Ionospheric  Time  Delay/Range  Error 

The  index  of  refraction  of  the  ionosphere  can  be  defined  by 

(1)  „  .  (x  .  13LS1)J 

when  electron  collision  frequency  and  magnetic-field  effects  are  ignored  ([U], 
p.  390).  Rote  that  the  refractive  index  is  directly  proportional  to  the  elec¬ 
tron  density  and  Inversely  proportional  to  the  square  of  the  frequency  of  the 
propagating  signal,  litis  index  of  refraction  causes  the  velocity  of  propaga¬ 
tion  to  be  less  than  that  in  free  space,  and,  hence,  Introduces  an  excess  time 
delay  or  range  error.  The  range  error,  in  meters,  can  be  expressed  by 
a 

(2)  R  -  V  N  ds 

***  o 

where  the  frequency,  f,  is  in  Hs  and  the  Integrated  electron  content  along  the 
ray  path  ia  in  electrons/**  [16).  Note  that  the  range  error  ia  inveraely  pro¬ 
portional  to  the  a qua re  of  the  frequency  of  the  propagating  signal  and  directly 
proportional  to  the  total  number  of  electrons  along  the  ray  path.  Typical 
range  errors  for  100,  200,  and  500  aHt  signals  are  illustrated  in  Figure  1. 
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Figure  1.  Range  Error  Versus  Total  Electron 
Content . 


SECTION  B  •—  THE  IONOSPHERE 


Production/Loss  Mechanisms 

The  earth's  ionosphere  is  produced  primarily  by  the  ionisation  of  neutral 
atmospheric  constituents  by  solar  radiation;  ionising  radiation  includes  par- 
ticles,  ultraviolet,  and  x-ray.  Because  this  report  is  concerned  mainly  with 
the  upper  regions  of  the  ionosphere  (the  R2  layer  and  the  "topside"),  it  is 
most  concerned  with  the  principal  source  of  ionisation  in  the  P  region:  solar 
Extreme  Ultraviolet  (EUV)  radiation  from  200-800  i. 

Electrons  aro  lost  by  recombination  with  positive  ions  or  by  attachment  to 
neutral  atoms. 

The  balance  between  production  and  loss  processes  in  the  F  region  is 
a  complex  function  of  variable  radiation  and  a  hlghly-variable  chemistry  of  the 
atmospheric  constituents  (see  Craig  (4),  Chapter  9;  Kelso  (14),  Chapter  3) • 

Transport  of  ionisation  along  magnetic  field  lines  is  also  an  important 
consideration  in  the  P  region.  Advectlon  of  electrons  into  or  out  of  an  area 
of  Interest  must  be  thought  of  ns  an  additional  gain/loss  mechanism;  this 
mechanism  is  thought  to  be  responsible  for  many  of  the  observed  changes  in 
electron  content. 

TEC  Measurements 

Figure  2  shows  a  typical  vertical  profile  of  daytime  electron  density  in 
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Figure  2.  A  Typical  Daytime  Ionospheric 
Profile,  Electron  ’'•ensity 
Versus  Height. 


the  ionosphere.  The  indicated 
D,  E,  and  P  layers  in  the  lower 
ionosphere  are  of  significance 
primarily  to  communicators  who 
rely  on  the  reflective  character¬ 
istics  of  the  ionosphere  in  the 
HP  (3-30  mHz)  region  of  the  spec¬ 
trum.  The  vertical  profile  must 
also  he  known  by  operators  of 
range -measuring  systems  with  a 
satellite  terminal  (or  target) 
actually  in  the  ionosphere. 
Systems  to  which  this  report  re¬ 
lates  are  those  using  VHF-UHF 
signals  which  pass  entirely 
through  the  ionosphere;  for  this 
reason,  the  entire  integrated 
electron  content  in  a  one-meter- 
square  vertical  column  from  the 
surface  to  an  altitude  of  a  few 


thousand  kilometers  will  be  treated  as  a  distinct  parameter,  analogous  to  a 


standard  meteorological  variable:  surface  pressure. 


Note  that  standard  ground-based  measurements  of  ionospheric  parameters  are 
limited  to  that  part  of  the  ionosphere  below  the  F2  peak  in  electron  density. 
Roughly  two-thirds  of  the  electrons  in  the  ionosphere  are  found  above  the  F2 
peak  and  hence  are  screened  from  conventional  ground-based  observations.  An 
important  relationship  exists  between  the  peak  electron  density,  NmP2,  and  the 
?'2-region  critical  frequency,  fQF2  (14),  p.  [%)x 

(3)  H/2  -  (££}  f*  •  1,24  *  10"’  (f0P2)# 

There  are  a  number  of  techniques  that  can  be  used  to  make  direct  observa¬ 
tions  of  total  electron  content.  Tie  great  majority  of  data  available  today 
have  been  obtained  by  analysing  the  rotation  of  the  place  of  polarisation  of 
VHP  signals  transmitted  by  satellite;  this  Faraday  rotation  technique  is  well 
documented  (8)  (23)  and  is  accepted  as  an  accurate  and  useful  measure  of  TEC 
up  to  an  altitude  of  a  few  thousand  kilometers.  Faraday  rotation  is  propor¬ 
tional  to  the  product  of  electron  density  and  the  magnetic  field  strength 
along  the  ray  path.  A  great  deal  of  literature  is  available  on  the  major  lim¬ 
itation  to  this  technique;  the  need  to  Include  effects  of  the  geomagnetic 
field  weighted  by  the  electron  density  profile  (1)  (26).  In  addition  to 
errors  Introduced  by  an  assumption  of  magnetic  field  strength  at  the  centroid 
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of  the  vertical  profile  along  any  penetrating  ray,  the  method  does  not  count 
electrons  at  higher  altitudes  where  the  strength  of  the  earth's  field  is  too 
low  to  contribute  to  the  rotation  effect.  The  data  used  in  the  development  of 
this  TEC  model  were  obtained  at  the  Air  Force  Cambridge  Research  Laboratories' 
Sagamore  Hill  (Mass.)  Radio  Observatory  by  the  Faraday  rotation  technique,  and 
the  model  is  intended  to  specify  and  predict  values  of  TEC  as  measured  through 
Faraday  rotation.  The  true  total  content  between  the  earth  and,  say,  a  geo¬ 
stationary  satellite,  will  be  perhaps  1056-1556  higher  than  the  values  given 
here. 


Vertical  Total  Electron  Content 

The  large-scale  morphology  of  mid-latitude  TEC  is  well  illustrated  in 
Figure  3  [16].  Marked  seasonal  and  diurnal  changes,  as  well  as  day-to-day 
fluctuations,  are  evident.  (For  a  description  of  the  complexities  of  worldwide 
TEC  morphology,  see  Klobuchar  [17].)  TEC  behavior  is  similar  in  many  respects 
to  fQF2,  the  critical  frequency  of  the  F2  layer;  this  fact  is  fully  exploited 
in  the  development  of  the  TEC  models  in  this  report.  The  qualitative  simi¬ 
larities  are  particularly  noticeable  in  the  seasonal  and  diurnal  variations; 
the  well-documented  (but  little  understood)  winter  or  seasonal  anomaly  in  fQP2 
also  shows  up  in  the  TEC  observations.  The  winter  anomaly  is  the  phenomenon 
of  higher  Northern  Hemisphere  mid-latitude  fQF2s  in  December  than  in  June 
(see  Rishbeth  and  Garriott  (21J,  Chapter  V). 

In  following  sections  of  this  report,  it  will  be  convenient  not  to  have  to 
work  with  an  ionosphere  such  as  the  one  shown  in  Figure  2;  an  ionosphere  hav¬ 
ing  a  single  electron  density,  constant  with  height,  is  much  easier  to  manipu¬ 
late.  An  artificial  Ionosphere  having  a  constant  electron  density  chosen 
equal  to  maximum  density  in  the  real-world  ionosphere,  NaF2,  ir  shown  in  Fig¬ 
ure  4.  The  thickness,  t,  of  such  a  rectangular  ionosphere  is  chosen  so  that 


(*> 


TEC 


where  N  is  the  single-valued  electron  density  and  TEC  is  the  vertical  total 
content  of  that  same  real-world  ionosphere. 

The  parameter,  t,  is  called  the  slab  thickness,  and  it  plays  a  major  role 
in  TEC  modeling  for  the  following  reason.  Note  that,  if  Equation  (3)  Is  sub¬ 
stituted  In  Equation  (4), 

TEC 

1.24  x  10“*  (f0F2)# 
or 

(5)  TEC  «  1.24  x  10**  (f0P2)“  t 
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TEC  can  thup  be  predicted  (or  speci¬ 
fied)  if  one  has  information  on  the 
behavior  f  the  _E2  critical  fre¬ 
quency  and  the  slab  thickness.  The 
climatology  of  fgF2  is  well-known; 
numerical  map  functions  for  the 
worldwide  variation  of  fifthly  me¬ 
dian  fQF2  are  prepared  routinely  by 
the  Environmental  Science  Services 
Administration  (ESSA),  and  both 
real-time  observations  and  short- 
range  forecasts  of  fQF2  are  avail¬ 
able  for  USAF  operational  applica¬ 
tions  from  the  Air  Weather  Service 
Aerospace  Environmental  Support 
Center.  The  ESSA  monthly  median 
predictions  [13!  are  used  in  the 
following  TEC  mcu.el  as  valid  cli¬ 
matological  descriptions  of  the 
lower  F  region  (fQF2). 

Slant  Total  Electron  Content 

Operational,  applications  of  TEC  models  will  seldom  involve  vertical  ray 
paths.  Because,  vertical  content  is  modeled  and  slant  extent  must  be  applied 
operationally,  correction  factor  must  be  applied  to  modeled  TEC  values. 

Taking  the  rectangular  ionosphere  as  an  example  (Figure  4),  it  can  be  seen  that 
the  slant  TEC  is  related  to  the  vertical  TEC  by 

<TECWt  -  (“‘'Ivertid  Ssc  8 

where  0  is  the  angle  between  the  ray  direction  and  the  zenith.  In  the  real- 
world  lonosphero  a  mean  value  of  Sec  9,  weighted  by  the  electron  density  pro¬ 
file,  must  be  determined.  This  weighted  mean  Sec  9  if  •normally  selected  at  an 
altitude  of  350  ’KJO  km.  The  zenith  angle  of  the  ray  at  the  ionospheric  pene- 
tration  point  (a t,  soy,  400  km)  is  a  function  of  the  elevation  angle  of  the 
ray  at  the  surface  [ 5 ] s 

Sin  9  «  Cos  e 
or 

(6)  Sec  9  -  [l  -  (syffiffirog)  cos*  e] 
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Figure  4.  An  Artificial  Ionosphere: 

Electron  Density  Versus 
Height . 
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Figure  5.  Geometry  of  ionospheric  Penetra¬ 
tion. 

Figure  5  illustrates  the  effect  of  elevation  angle  on  penetration-point  senith 
angle j  note  the  maximum  8  at  an  elevation  angle  of  0°.  8  decreases  to  0e  at 

an  elevation  angle  of  90° . 

Values  of  Sec  8  for  varying  elevation  angles  can  be  determined  from  Fig¬ 
ure  6.  Note  the  maximum  value  of  Sec  8  *  2.95  for  a  0° elevation  angle.  Thus, 
the  slant  TEC  can  be  nearly  three  times  greater  than  the  vertical  TEC  at  the 
ionospheric  penetration  point. 


\ 
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Figure  6.  Sec  9  Versus  Elevation  Angle. 


SECTION  C  —  SOLAR  RADIATION 

The  EUV  radiation  that  produces  the  upper  ionosphere  is  emitted  from  the 
hot,  tenuous  gaBes  in  the  lower  solar  corona.  Direct  observations  of  solar  EUV 
can  be  made  only  from  satellite  platforms  because  the  radiation  is  completely 
absorbed  by  the  atmosphere  before  it  can  reach  the  earth's  surface. 

Solar  radiation  In  radio  wavelengths  is  emitted  by  the  same  region  of  the 
sun's  atmosphere,  and  observations  of  solar  radio  noise  have  been  used  to  infer 
changes  in  EUV  radiation  (COESA  (3],  Part  3)>  Solar  radio  noise  at  a  wave¬ 
length  of  10  7  cm  has  been  observed  continuously  since  the  late  1940s;  this 
measure  of  the  sun's  radio  output  (in  watts/ma/Hz)  has  been  successfully  used 
in  models  of  neutral  density  at  satellite  altitudes  as  a  measure  of  EUV  heat¬ 
ing  of  the  upper  atmosphere.  The  radio-flux  parameter,  F10  (units  of  10"aa 
watts/m*/Hz),  is  observed  by  the  National  Research  Council  of  Canada,  Ottawa, 
and  by  the  Air  Forte  Cambridge  Research  Laboratories,  Sagamore  Hill,  Massa¬ 
chusetts.  Data  are  readily  available  to  operationkl  and  research  agencies 
from  the  Air  Weather  Service  Space  Environmental  Support  System.  F10  is  used 
in  this  study  to  model  solar  EUV  radiation.  References  to  F10  should  also  be 
interpreted  as  references  to  solar  EUV. 

It  is.  possible  to  consider  the  daily  values  of  F10  as  consisting  of  two 
components  (22] :  (1)  a  basic  component,  associated  with  a  hypothetical  sun- 

spot-free  sun,  and  (?)  an  active-region  component  (also  known  as  the  plage 
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component  or  the  slowly-varying  component),  associated  with  the  higher  coronal 
temperatures  and  densities  above  individual  centers  of  activity  on  the  sun's 
surface.  The  basic  component  varies  with  the  11-year  solar  cycle,  and  the 
active-region  component  shows  both  strong  day-to-day  variations  and  a  27-day 
periodicity  corresponding  to  the  rotation  of  active  regions  across  the  visible 
disk  of  the  sun.  Typical  monthly  mean  values  of  P10  for  the  last  20  years  are 
shown  in  Figure  7.  A  smoothed  estimate  [10]  of  the  background  component,  B, 
obtained  during  Sunspot  Cycle  number  19  (195^-1964)  is  also  shown  in  Figure  7* 
The  background  component  varies  from  68-70  units  at  solar  minimum  to  an  esti¬ 
mated  130  units  at  solar  maximum.  The  value  of  B  at  solar  maximum  may,  indeed, 
vary  from  sunspot  cycle  to  sunspot  cycle 5  however,  for  the  purpose  of  this 
study  it  is  assumed  that  the  solar  radio  (EUV)  flux  may  be  described  by  two 
components:  (1)  a  background  value,  B  (nearly  constant  during  colar  maximum 
or  minimum  years),  and  (2)  an  active  region  contribution,  F10  -  B,  which  varies 
with  the  movement,  growth,  and  decay  of  centers  of  activity.  Thus: 

FLUX  -  B  +  (F10  -  3) 


SECTION  D  —  TEC  MODEL  I 


Assumption 

The  basic  premise  of  the  TEC  model  is  that  peak  dally  TEC  values  at  a 
"typical"  Northern  Hemisphere  mid-latitude  location  can  be  related  to  a  measure 
of  the  solar  EUV  radiation;  TEC  changes  associated  with  geomagnetic  disturb¬ 
ances  are  not  considered.  For  this  effort,  it  is  assumed  that  the  peak  daily 
TEC  value  for  a  "standard"  mid-latitude  station  can  be  specified  by: 

(7)  TECp  «  CQ  +  Cx  (no  -  B) 

TEC  Components 

The  first  term  of  Equation  (7),  CQ»  represents  the  TEC  contribution  of  the 
solar  background  radiation,  B.  This  term  includes  the  seasonal  changes  in  TEC 
caused  by  primarily  terrestrial  ratner  than  solar  sources.  is  illustrated 
in  Figure  8.  This  particular  curve  was  derived  from  observed  TEC  values  dur¬ 
ing  sunspot  maximum;  suggest Iona  for  similar  functions  applicable  to  other 
stages  of  the  tolar  cycle  are  given  in  Section  0. 

The  second  term  in  Equation  (7)  represents  the  TEC  contribution  of  solar 
active-region  radiation.  A  number  of  researchers  (2)  (7)  [12]  have  reported 
both  linear  and  logarithmic  relationships  between  TSC  and  either  solar  radio 
flux  or  smoothed  sunspot  number;  *11  report  eeasonel  variations  in  the 
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Figure  8.  TEC  Basic  Component,  CQ. 


constants  of  proportionality.  Because  seasonal  effects  are  modeled  in  the 
basic  component  of  TEC,  the  relationship  between  TEC  and  solar  active-region 
radiation  is  assumed  to  be  constant  throughout  the  year.  {A  test  of  TEC  Model 
I,  albeit  on  data  from  which  the  model  was  derived,  snowed  no  improvement  in 
TEC  specification  when  the  value  of  C^  was  adjusted  seasonally.)  A  represent¬ 
ative  value  of  the  constant  of  proportionality,  C^,  was  selected  from  an 
analysis  of  the  relationships  reported  in  the  literature;  a  value  of  0.03  was 
chosen  for  use.  A  further  refinement  of  the  second  term  in  Equation  (7)  is 
necessary  to  remove  the  "noisiness"  of  day-to-day  fluctuations  in  PlOj  there¬ 
fore,  P10  is  smoothed  by  the  use  of  a  mean  value.  For  the  specification  of  a 
monthly  mean  peak  TEC  value,  the  observed  or  predicted  monthly  mean  P10  is 
used.  For  the  specification  or  prediction  of  daily  peak  TEC  values,  a  5-day 
running  mean  of  F10  (TOJ^)  is  used:  computed  for  days  01  through  05  is 

used  to  specify  the  peak  TEC  values  for  05  or  to  predict  the  peak  value  for 
06.  The  daily  peak  TEC  value  is  giv  n  by 

TECp  -  CQ  +  0.03  -  B) 

where,  during  sunspot  maximum,  CQ  is  determined  from  Plgujre  8  and  B  is  130 
flux  units,  giving 

(8)  TECp  -  C0  +  0.03  (TOJg  -  130) 

Diurnal  Variation  of  TEC 

A  diurnal  curve  it  instructed  around  the  specified  or  predicted  value  of 
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TECp.  The  diurnal  variation  for  the  standard  mid-latitude  location  is  empiri¬ 
cal  in  nature  and  is  based  on  observations  provided  by  AFCRLj  the  "standard" 
mid-latitude  station  is  chosen  to  duplicate  the  AFCRL  Sagamore  Hill  Radio  Ob¬ 
servatory,  Massachusetts,  with  an  ionospheric  penetration  point  near  4o°N  75°W. 
Five  key  values  of  TEC  versus  local  time  are  selected  and  are  merely  connected 
with  straight  lines.  The  seasonal  variation  in  the  shape  of  the  curve  is  pro¬ 
duced  by  shifting  the  local  time  of  the  peak  value.  The  construction  of  the 
curve  is  described  below,  and  is  illustrated  in  Figure  9. 


Fleur*  9.  Seasoned  Variation  of  Diurnal  Curve  Construction. 
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TEC  MOD  I 

Predict  diurnal  TEC  curve  for  standard  mid-latitude  site  (40°N  ?5°W)  for 
day  d6: 


a.  Predict  peak  TEC,  TECpi 

(1)  Determine  basic  component  of  TECp,  C0j  use  Figure  8. 

(2)  Determine  5-day  running  mean  of  F10,  FlO^,  for  days  D1  through  D5. 

(3)  Using  Equation  (8),  TECp  **  CQ  +  0.03  (FlS^  -  130). 

b.  Specify  local  time  of  TEC  : 


1400  LST 
1600  LST 
1800  LST 


Winter  (Nov-Feb) 

Equinox  (Mar-Apr,  Sep-Oct) 

Summer  (May -Aug) 

c.  Specify  TEC  *  TEC  -  0.5  x  1017  at  1100  LST  and  (Peak  time  +  2  hrs). 

y 

d.  Specify  TEC  at  2100  and  2400  LST:  See  Figure  10  (Note:  x  is  the  peak 
solar  zenith  angle  at  4o°N  on  D6.) 


Figure  10.  Factor*  for  TSC  Specification  at 
2100,  2400  LST. 
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(1)  TEC2100  =  TECp  (.64  Cos1-8  x).  Use  Figure  10a. 

(2)  TEC2i;00  ■  TECp  (.45  Cos1-6  x).  Use  Figure  10b. 

e.  Specify  TEC  at  sunrise  =  0.5  x  1017 

Winter  0600  LST 
Equinox  0500  LST 

Summer  0400  LST  „  , 

TABLE  1 

Example  of  TEC  MOD  I  j  DATA 


A  specification  of  the  monthly  median 
diurnal  curve  for  4o°N  75°W,  April  1969*  is 
illustrated  below.  Using  the  data  in  Table 
1,  the  observed  monthly  mean  of  F10  is  used 
in  place  of  FlO^  which  would  be  used  for 
daily  TEC  predictions. 

Step  1.  CQ  -  2.8  x  1017 

Fl530  -  155 

TEC  -  [2.8  +  0.03  (25) J  x  10* * 

-  3.6  x  101V 

Step  2.  Tine  of  TECp  »  1600  LST 

Step  3.  TEC  at  1100,  1800  LST  »  3.1  *  10* 7 

Step  4.  TEC  at 

2100  *  3.6  *  10*7  (.501 

-  1.8  x  10*7 

2400  -  3.6  x  10*7  (.35) 

»  1.4  x  10*’ 

Step  5.  TEC  at  sunrise  (0500  LST) 

-  0.5  x  10*7 

Step  6.  Points  aro  connected  with  straight 
lines;  a  comparison  with  observed 
data  is  shown  in  Figure  11. 


DATE 

(1969) 

F10 

m5 

Co 

TECp 

27  Mar 

178 

.  « 

28 

178 

— 

•  • 

-- 

29 

182 

... 

.. 

30 

18.3 

•  • 

31 

186 

181 

-- 

-- 

1  Apr 

189 

184 

3.1 

4.6 

2 

191 

186 

3.1 

4.7 

3 

190 

188 

3.0 

4.7 

4 

177 

187 

3.0 

4.7 

5 

176 

185 

3.0 

4.7 

6 

163 

179 

3.0 

4.7 

7 

I55 

172 

3.0 

4-5 

8 

i58 

164 

3.0 

4.3 

9 

144 

157 

2.9 

3-9 

10 

149 

152 

2.9 

3.7 

11 

151 

149 

2.9 

3.6 

12 

156 

150 

2.9 

3.5 

13 

173 

155 

2.9 

3.5 

i4 

179 

162 

2.8 

3.6 

15 

181 

168 

2.8 

3.8 

16 

167 

171 

2.8 

4.0 

17 

155 

171 

2.8 

4.0 

18 

147 

166 

2.8 

4.0 

19 

146 

159 

2.8 

3.9 

20 

149 

153 

2.7 

3.6 

21 

157 

151 

2.7 

3.4 

22 

148 

149 

2.7 

3.3 

23 

144 

149 

2.7 

3.3 

24 

145 

149 

2.7 

3.3 

25  , 

148 

148 

2.7 

3.3 

26 

144 

146 

2.6 

3.1 

27 

135 

143 

2.6 

3.1 

28 

141 

2.6 

3.0 

29 

128 

138 

2.6 

2.9 

30 

129 

2.6 

2.8 

15 
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Figure  11.  Specified  end  Observed  TEC, 
April  1969. 

SECTION  E  —  TEC  MODEL  II 


Assumption 

TSC  MOD  II  is  designed  to  generate  »  diurnal  curve  for  locations  other 
than  the  "standard"  aid-latitude  station.  The  assumption  necessary  to  convert 
the  standard  curve,  TECg(t),  to  any  other  location,  TEC^(t),  is  that  the  slab 
thickness  at  any  given  local  tine  is  the  sane  the  world  overt  ■  t#.  The 
few  available  comparisons  of  slab  thickness  at  different  stations  U5)  [24} 
indicate  that  the  assumption  is  reasonable  over  middle  latitudes.  The  assump¬ 
tion  is  not  valid  at  equatorial  latitudes,  particularly  during  sunlit  hours. 

Derivation 

Applying  the  assumption  of  identical  diurnal  slab-thickness  variations  at 
all  points, 

yt)  ■  y*)* 

By  Equation  (4),  then,  at  any  given  local  time, 
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N. 

(9)  TECi  =  TECg  TJi 
Substituting  for  N  from  Equation  (3) 

(V2i>a 


TEC1  =  TECS 


(VV 


or 


(10)  TECj,  *  TECS  •  L 


where 

(U>  L  .  ($£)’ 

o  s 


Application 

The  conversion  factor,  L,  can  be  computed  from  an  analysis  of  median  fQP2 
predictions.  One  can  then  specify  TOC  at  any  point  at  a  «iven  local  time  by 
applying  the  L-f actor  to  the  standard  TEC  prediction  at  that  local  time.  Ta¬ 
ble  2  shows  an  example  of  L- factors  computed  for  longitudes  15°E-l65°E  for  the 
month  of  April  1969,  each  L- factor  is  valid  at  1100  LST  at  the  different  long¬ 
itudes. 


TABLE  2 

L-Pactors  for  1100  LST,  April  1969. 


East  Longitude 


North 

Latitude 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 

Mean 

70 

0.7 

0.7 

0.8 

0.9 

0.9 

0.8 

0.8 

0.7 

0.6 

0.6 

0.5 

0.7 

60 

0.6 

0.9 

l.o 

1.2 

1.2 

1.2 

1.3 

1.0 

0.9 

0.8 

0.8 

1.0 

50 

1.1 

1.1 

1.2 

1.4 

1.5 

1.5 

1.5 

1.4 

1.3 

1.2 

1.2 

1.3 

40 

1.5 

1.5 

1.6 

1.7 

1.8 

1.9 

1.9 

1.8 

1.8 

1.8 

1.7 

1.7 

30 

1.9 

2.2 

2.2 

2.3 

2.3 

2.3 

2.3 

2.4 

2.3 

2.4 

2.2 

2.3 

20 

2*3 

2.3 

2.5 

2.8 

2.6 

2.4 

2.5 

2.7 

2.5 

2.4 

2.3 

2.5 
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The  TEC  at  1100  LST  at  each  point  in  the  table  would  be  computed  by  multi¬ 
plying  the  standard  TEC  from  Figure  10  at  1100  LST  (3.1  *  1017  e/ms )  by  the 
appropriate  L-factor,  as 

TECi,  1100  LST  =  TECs,  1100  LST  *  Li,  1100  LST 
It  is  sometimes  convenient  to  require  L  to  be  of  the  form 

L  =  A  exp  (-bcp) 

where  cp  is  geographic  latitude.  As  an  example  of  this  form  of  L,  the  latitu¬ 
dinal  variation  of  L  at  a  given  local  time  and  given  longitude  (or  range  of 
longitudes)  is  plotted  and  a  best-fit  exponential  for  the  406N-60°N  region  is 
determined.  Figure  12  shows  the  mean  Eastern  Hemisphere  L-factor  versus  lati¬ 
tude  from  Table  2j  the  best-fit  exponential 

L  =  5.1  exp  (-.027  <P) 

is  also  shown.  The  illustrated  curve  is  an  exceptionally  good  fit,  and  is  not 
expected  in  all  cases,  particularly  at  low  latitudes. 


Figure  12.  L- Factor  Versus  Latitude,  1100  LST. 
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Example  of  TEC  MOD  II 

The  specified  monthly  median  TEC  diurnal  variation  for  Hawaii  is  derived 
for  March  1967. 


Step  1.  The  d’urnal  curve  for  the 
standard  station  is  gener¬ 
ated  as  in  MOD  I,  using 
™30  *  161.  The  standard 
diurnal  curve  is  shown  in 
Figure  13c. 

Step  2.  Diurnal  f  F2  curves  for  both 
Hawaii  and  the  standard  sta¬ 
tion  are  shown  in  Figure  13a* 

Step  3.  The  conversion  factor,  L, 

determined  by  tne  square  of 
the  ratio  of  the  fQf2  curves, 
is  ehown  In  Figure  l?b. 

a  ft*  *  P 

L'  'W 

Step  <».  TEC  *  7ECe  •  L 


Figure  13*  Example  of  ?SC  HOD  11. 

The  Modeled  curve,  defined  by  the  product  of  the  standard  TEC  curve  in 
Figure  13a  and  the  L-curve  In  Figure  13b,  is  ehown  in  Figure  13c. 


Technical  Report  234 


June  1970 


t 

SECTION  P  —  TEC  MODEL  III 


A  prime  deficiency  of  MOD  II  is  the  artificial  control  of  the  worldwide  di¬ 
urnal  curves  by  the  empirical  constraints  placed  on  the  standard  diurnal  curve. 
A  more  general  form  of  diurnal  variation  is  obtained  by  allowing  each  loca¬ 
tion's  TEC  to  vary  as  the  local  (f0F2)a.  MOD  II  assumed  equivalent  slab  thick¬ 
nesses  at  all  points  at  the  same  local  time;  an  unspecified  diurnal  variation 
in  t  was  allowed.  MOD  III  must  as¬ 


sume  an  absolutely  constant  slab 
thickness.  Because  t  is  known  to 
have  a  considerable  diurnal  variation, 
this  assumption  Introduces  rather  ob¬ 
vious  errors.  Slab  thickness  is 
chosen  to  be  the  standard-station 
slab  thickness  at  the  time  of  peak 


Derivation 


In  this  case,  the  peak  TEC  at  lo¬ 
cation,  1,  is  predicted  similarly  to 
MOD  II,  but  without  the  artificial 
specification  of  the  time  of  peak 
TEC.  Each  station's  diurnal  TEC 
curve  is  normalised  to  its  (fQP2)a 
diurnal  peak  value.  Assume 

Ti  “  Tsp  " 

Prom  Equation  (4), 


(12) 


TECl 

TEC1 


or 


(13)  TECj  -  TECfp  e  J 


Figure  14.  Example  of  TEC  MOD  HI. 
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To  predict  diurnal  TBS  curve  for  any  location,  i,  for  day  D6: 

a.  Predict  standard  ^eak  TEC,  TECsp: 

(1)  Determine  bas  •  component  of  TEC  ,  C  ;  use  Figure  8. 

“Jr 

(2)  Determine  5-$gr  running  mean  of  F10,  f16^,  for  days  D1  through  D5. 

(3)  TECsp  =  CQ  +  0,03  (FTO5  -  130). 

b.  Determine  factor,  J„  relating  TEC  to  TEC.: 

(1)  Determine  peak,  value  of  fQF2  diurnal  variation  at  standard  loca¬ 
tion:  40° N  75° W. 

(2)  Determine  diurimi  variation  of  f0F2  at  location,  1. 

,f  F2.  ( t]k2 

(3)  J(t)  =  (~f ) 

o  sp 

c.  TECi(t)  -  TDCsp  °  J(t| 

TEC  MOD  III  gives  a  reasonable  estimate  of  mid-latitude  diurnal  TEC  varia¬ 
tions  under  geomagnetically  %^t  conditions.  Derivations  from  absolute 
(measured)  values  are  large  aS  low  latitudes  at  all  times  of  day,  and  are 
largest  at  mid-latitudes  during  the  night  and  post  sunrise  hours.  The  shapes 
of  the  diurnal  curves  are,  however,  more  realistic  at  all  latitudes  than  those 
in  MOD  II.  While  TEC  absolute,  value  errors  are  greater  in  MOD  III,  longitu¬ 
dinal  TEC  gradients  inferred  from  the  diurnal  curve  are  more  accurate.  Im¬ 
provements  *n  this  model  must  await  additional  knowledge  of  the  worldwide  mor¬ 
phology  of  slab  thickness.  The  Appendix  ’'tains  a  general-purpose  computer 
program  for  TEC  MOD  III  predictions  based  on  ESSA/ITS  ionospheric  predictions. 


Example  of  TEC  MOD  III 

The  specified  monthly  median  TEC  diurnal  variation  for  Hawaii  is  again  de¬ 
rived  for  March  1967. 

Step  1.  The  standard  peak  TEC  value  is  determined  as  in  MOD  I. 

TEC  -  CQ  +  0.03  (m#  -  130)  »  3.2  +  0.03  (31)  «  4.1  x  10i?  e/«a 

Stop  2.  Determine  factor  J: 

/f  F2,  Na 

o  sp 

The  predicted  monthly  media  f.F2  curve  for  Hawaii  and  the  maximum  f_?2  for 
the  standard  location  are  shown  in  Figure  l4a.  J  is  shown  in  Figure  l4b, 


Step  3.  TECi  -  TECap  «  J 

The  product  of  J  (Figure  14k)  and  the  standard  peak  TEC  is  shown  in  r  mu 
14c .  '  ■> 
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SECTION  G  —  MODIFICATIONS 


Geomagnetic  Disturbances 

The  behavior  of  mid-latitude  TEC  during  geomagnetic  disturbances  has  so  far 
resisted  all  attempts  at  modeling.  A  few  general  rules  have  been  extracted 
from  TEC  observations  [11]  [20]: 

a.  TEC  increases  above  normal  quiet -day  values  during  the  first  day  of  a 
sudden-commencement  (SC)  disturbance;  the  increase  is  typically  greatest  dur¬ 
ing  pre-sunset  hours  on  the  first  day  of  the  disturbance.  Disturbed  values  of 
TEC  may  be  double  the  typical  quiet-day  values. 

b.  TEC  is  generally  below  quiet-day  values  for  the  remainder  of  the  dis¬ 
turbance;  the  largest  negative  percentage  deviations  are  observed  between 
0300-0600  LST.  Mendillo,  in  a  detailed  study  of  22  geomagnetic  disturbances, 
reports  apparent  seasonal  as  well  as  storm-time  characteristic  TEC  behavior. 

c.  There  are  no  consistent  relationships  observed  between  standard  geo¬ 
magnetic  indices  (Ap,  K)  and  TEC  variations;  TEC  is  apparently  better  related 
to  storm  phase. 

d.  There  are  indications  that  f  ?2  changes  precede  TEC  changes  by  a  few 
hours . 

Slab  Thickness 

It  is  clear  that  one  need  go  no  further  than  Equation  (5>)  to  model  TEC  if 
the  slab  thickness  is  a  known  quantity.  Klobuchar  [IS]  is  developing  an  ex¬ 
pression  for  mid-latitude  daytime  t  (at  Sagamore  Hill,  Massachusetts)  of  the 

form 

(14 )  t  ®  a  +  B  Sin  ^  (h  -  C)  +  E  Sin  (i  -  F) 

;uhere  A,  R,  C,  E,  and  P  are  constants,  h  is  the  local  time,  and  d  is  day  of 
\r.  If  this  expression  can  be  extended  to  other  latitudes  (and  longitudes), 
hen  reel-time  observations  (or  she "t- range  predictions)  of  f£)‘,'2  could  be  used 
specify  (or  predict)  TEC  with  an  expectation  of  high  accuracy. 

v  Until  further  definition  of  Equation  (i4)  can  be  attained,  it  may  be  pos¬ 
sible  to  improve  the  applicability  of  TEC  MOO  IX  or  MOO  III  to  low  latitudes 
by  modifying  the  assumption  of  c.  urt  slab-thickness.  Observations  of 
monthly-median  *lab  thickness  at  Sagamore  Hill  suggoat  a  relationship  between 
midday  *  and  the  noon-ti»a  solar  nenitU  angle  of  the  form 

t  «  A  +  B  Cos*  y 
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The  Cos0  x  relationship  looks  particularly  good  during  the  summer  months  and  is 
illustrated  in  Figure  15 .  (See  Goodman  [9]  for  a  presentation  of  a  similar  re¬ 
lationship  between  t  and  Cos  x  which 
is  extended  to  an  expression  for  TEC 
as  a  function  of  Cos  X. )  We  may  as¬ 
sume  that  similar  relationships  exist 
between  slab  thickness  and  x  at  other 
latitudes;  relationships  may  then  be 
developed  among  slab  thicknesses  at 
different  locations  using  Cos  x  or 
Cos2  x.  Such  relationships  could 
then  be  applied  to  Equation  (14)  to 
specify  r  at  any  point. 

Solar  Cycle  Variation  of  C0 

The  term,  CQ,  in  Equation  (7)  de¬ 
scribes  that  component  of  TEC  related 

r 

to  the  meteorology  of  the  earth' s  upper 
atmosphere  and  to  the  background  com¬ 
ponent  of  solar  ultraviolet  radiation. 

CQ  should,  then,  vary  with  the  solar 
cycle.  By  combining  data  reported  by 
Youakim  and  Rao  [24]  with  Sagamore 
Hill  observations,  a  relationship  be¬ 
tween  the  minimum  annual  (summer)  C  Figure  15 .  Midday  Slab  Thickness 

and  solar  background  radiation,  B,  is  £sin!  ^t^SnloUr 

estimated  to  be  Zenith  Angle. 


(15)  MIN  CQ  =  0.13B 

The  amplitude  of  the  semiannual  CQ  variation  is  estimated  to  be 

(16)  MAX  CQ  -  MIN  CQ  =  0.115B 

Using  Equations  (15)  and  (16),  it  is  possible  to  estimate  the  variation  irt, C0 
over  the  solar  cycle.  Suggested  CQ  curves  for  different  phases  of  the  solar 
cycle  are  presented  in  Figure  16.  The  shape  of  these^ curves  may  also  vary 
with  solar  cycle.  One  particular  aspect  of  the  suggested  CQ  variation  should 
be  watched  closely  by  operational  agencies:  the  change  from  characteristic 
equinoctial  to  summer  behavior  is  often  observed  to  occur  over  a  short-time 
span  (say,  on  the  order  of  weeks).  The  gradual  seasonal  changes  in  C0  indi¬ 
cated  in  Figures  8  and  16  may  have  to  be  adjusted  to  reflect  a  more  rapid 
changeover  from  one  season  to  another. 
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Figure  16.  Solar  Cycle  Variation  of  CQ . 


SECTION  H  —  CONCLUSIONS 


The  TEC  models  presented  in  this  report  were  generated  in  recognition  of 
an  immediate  need  for  operationally-oriented  data.  The  use  of  the  models  re¬ 
quires  a  minimum  of  resources:  all  are  suitable  for  hand  calculations,  and 
requisite  input  data  are  readily  available.  Key  points  are  reiterated: 

a.  Assumptions  made  in  the  derivation  of  the  models  have  introduced  sig¬ 
nificant  limitations. 

b.  The  models  attempt  to  duplicate  Faraday- rotation  data,  and  hence  ac¬ 
count  for  approximately  9056  of  the  true  TEC  between  the  earth  and  a  geosta¬ 
tionary  satellite. 

c.  The  models  are  applicable  to  an  undisturbed  Northern  Hemisphere 
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mid-latitude  ionosphere  near  the  maximum  of  the  solar  cycle.  They  are  not 
valid  near  the  equatorial  anomaly,  in  the  auroral  zone,  or  during  geomagnetic 
disturbances . 

The  accuracies  that  may  be  expected  are  as  follows: 

a.  Standard  station  monthly -median  peak  values  should  be  accurate  to  with¬ 
in  about  5^5 

b.  Standard  station  daily  peak  values  should  be  accurate  to  within  about 
15 

c.  Specification/prediction  of  peak  daily  values  at  an  arbitrary  location 
should  be  accurate  to  within  about  25$. 

This  report  is  concluded  with  an  opinion:  that  it  has  underscored  the  need 
for  far  more  investigation  than  has  been  accomplished  to  date.  An  operation¬ 
ally  meaningful  specification  of  a  significant  environmental  parameter,  such  as 
TEC,  must  be  valid  everywhere  the  Air  Force  has  to  function  —  and  this  cer¬ 
tainly  includes  auroral-zone  and  equatorial  latitudes.  It  appears  that  the 
future  of  operational  TEC  specification  lies  in  relating  observable  ionospher¬ 
ic  parameters,  such  as  fQF2,  to  functions  such  as  suggested  by  Klobuchar 
[Equation  (14)]  that  could  be  applied  anywhere  in  the  world.  All  of  which 
serves  to  emphasize  the  need  for  expanded  research  efforts. 
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APPENDIX 


The  computer  subroutines  listed  in  this  appendix  are  designed  to  calculate 
the  value  of  L  in  Equations  (10)  and  (11)  or  J  in  Equations  (12)  and  (13)  for 
use  in  the  calculation  of  total  electron  content.  This  program  relies  heavily 
upon  the  work  described  in  ESSA  Technical  Report  ERL107-ITS  75>  May  1969, 
Advances  in  Ionospheric  Mapping  by  Numerical  Methods.  Subroutines  GK,  SICOJT, 
DKSICO,  DKGK  are  extracted  directly  from  the  report;  the  only  changes  made 
were  those  necessary  to  convert  them  to  conventional  FORTRAN  IV.  The  subrou¬ 
tines  TEC  and  SAGMAX  will  be  described  below. 

Subroutine  TEC  will  compute  the  value  for  any  point  at  any  time.  The  sub¬ 
script  i  refers  to  the  particular  point  desired;  the  subscript  SAG  refers  to 
the  reference  value  at  Sagamore  Kill  Observatory.  The  subroutine  calling  se¬ 
quences  is  (K,  U,  ALAT,  ALONj  AMaGLT,  TIME,  RATIO).  The  parameters  are  de¬ 
fined  as  follows: 

K,  U  Parameters  contained  in  ESSA  fQ(F2)  coefficients.  This 

coefficient  deck  may  be  read  with  the  READU  subroutine 
described  in  ERL107-ITS  75- 


ALAT,  ALON  Latitude  and  longitude  of  the  desired  point  in  degrees. 

Latitude  is  positive  north,  longitude  positive  east. 

AMAGLT  The  modified  magnetic  dip  for  the  desired  location.  This 

parameter  may  either  be  calculated  as  described  in 
ERL107-ITS  75  or  extracted  from  Figure  2  of  that  report. 
AMAGLT  is  expressed  in  degrees  with  positive  being  north. 

TIME  Expressed  in  universal  hour  angle  from  -180  to  +180 

degrees. 


RATIO  Output  parameter, 

rfo(F2>i  ia 

* 

Subroutine  SAGMAX  computes  the  peak  daily  value  of  [fQ(F2)gAG]a .  This 
parameter  is  used  as  the  denominator  in  the  variable  RATIO,  output  from  TEC. 
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SUBROUTINE  TEC ( K,U,AL AT, ALCN.AMAGLT, TIME, RATIO! 

DIMENSION  K(10),U(17,76),C(ti)*S(ti),G(76),D(76) , COORD ( 3 ) 

COORD f I )  =  AMAGLT 

COORD ( 2 >  *  ALON 

COORCO)  =  ALAT 

NT  =  K<9)  ♦  1 

CALI  GK ( K , COCRC , G ) 

CALL  SIC0JT(8,C,S,TIME> 

CALL  DKSICC(NT,K(10),U,S,C,D) 

CALL  DKGK(NT,G,0,ANS) 

ANS  *  ANS«ANS 

CALL  SAGMAX ( K ,U, $MAX I 

RATIO  =  ANS/SMAX 

RETURN 

END 


subroutine  SAGMAX(K,U,SMAX ) 

DIMENSION  K(10),U(17,76)«C(8)*SI8),G(76),D(76), SAG ( 3  I 
DATA  SAG/55.0. -75.0, *0.0/ 

NT  -  KI9I  ♦  l 
CALL  GKIK.SAG.G) 

eto  *  0.0 

DU  10  1*1,24 
T  *  1 5* l ~IB0 
CALL  SIC0JTI8.C.S.T) 

CALL  DKSICU*NT,KC10),U,S.C,0> 

CALL  CKGK  ( NT  ,  C  ,1'  ■ ANS I 

ANS  *  ANS* ANS 

IE 1 1)10.1  T.  ANS I  BIG  •  ANS 

10  continue 

SMAX  »  BIG 
RETURN 

INC 


SUBROUTINE  OxSICCt  MX.LH.U.S.C.UI 
0 1  RE  NS  1  UN  Utl?,7M,Ct8)«StS),Dt?b) 

00  5  X*l,RX 
DUI  «  Uli.M 
00  5  l*l.LM 

S  DIM  .  tilKI  ♦  U t  ?*L  « K ) *S  I L  I  *  U I  /*L  *  l ,  X  I  »C  I L  I 

kE  turn 

tNi1 


SU8RLUTINE  UXCMMX,U,U,ANS» 
OlRENSiiJN  Gl  76). 01  76) 

ANS  *  C.O 
fcl*  *  X*  l  »RX 

5  ANS  *  ANS  »  OUl.tlXl 
RETURN 
END 
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SUBROUTINE  GKIK.C.G) 

COMPUTE  COORDINATE  FUNCTIONS,  611),  1*1, ....K+l 

C( 1 )=MQDIFIED  LATITUDE.  C I  2  ) ,  CI3>,  =GEOG.  LONGITUDE  AND  LATITUDE 
DIMENSION  KI1),  Cll),  Gil) 

OATA  DR/0. 017453293/, N/8/ 

X»0R*C 
Y*C I  2 ) *DR 
Z*DR*C<3) 

K0*K 

SX*SIN( X ) 

G(2 )*SX 
G»l. 

IF  IKQ  ,EQ •  1)  GG  TO  12 
4  DO  10  1*2, K0 
10  Gmi>*SX*G!I  ) 

12  KDIF«M2)-K0 

IF  (KOIF  ,NE.  0)  GO  TO  16 
80  RETURN 
lb  J*1 

cxi-cosm 

CX*C*l 

T»Y 

18  KC«Kljm 

GIKC-2 ) *CX*CGSI T  ) 

G(KC-l)«CX«SIhm 
20  KN*K I J, 1 ) 

IF  IKDIF  .EO.  2)  GO  TU  28 
IF  IXC.GT.  75)  GO  TO  28 
00  22  l*KC,RN,2 
Gm*iX«Cl  1-2) 

22  Gi l»l)«SX*GI 1-1) 

28  IF  U  .EG.  Nl  GO  TO  80 
30  KOIF  *KI J*2 ) -RN 

IF  I KOIF  .£0.  0)  GO  TO  80 
32  CX*CX*CX1 
J«4*l 
F  J*  J 
T  *F  J*V 
GO  Tf  18 
ENO 


SUBROUTINE  SltOJHl ,C, 4, A) 

COMMUTE  SIM  JT  I,  COSUTI,  ...  ,L  FOR  ANCLE  A 

DIMENSION  Ctl),  Sill 

T  *0.01 745 J29*A 

c«cusm 

S'SIMTI 

00  19  1*2,1 

cm*c«c<i-i>-s«si  i-n 
10  sm*c*stt-u*s*c(i-i) 

RETURN 
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